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Abstract: We demonstrate that tunable attractive (bonding) and repulsive 
(anti-bonding) forces can arise in highly asymmetric structures coupled 
to external radiation, a consequence of the bonding/anti-bonding level 
repulsion of guided-wave resonances that was first predicted in symmetric 
systems. Our focus is a geometry consisting of a photonic-crystal (holey) 
membrane suspended above an unpatterned layered substrate, supporting 
planar waveguide modes that can couple via the periodic modulation of the 
holey membrane. Asynmietric geometries have a clear advantage in ease 
of fabrication and experimental characterization compared to symmetric 
double-membrane structures. We show that the asymmetry can also lead 
to unusual behavior in the force magnitudes of a bonding/antibonding 
pair as the membrane separation changes, including nonmonotonic de- 
pendences on the separation. We propose a computational method that 
obtains the entire force spectrum via a single time-domain simulation, by 
Fourier-transforming the response to a short pulse and thereby obtaining the 
frequency-dependent stress tensor. We point out that by operating with two, 
instead of a single frequency, these evanescent forces can be exploited to 
tune the spring constant of the membrane without changing its equilibrium 
separation. 
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1. Introduction 

The coupling of the electromagnetic field and matter can give rise to optical forces on oth- 
erwise neutral objects and optomechanical coupling via radiation pressure or gradient 
forces in nanophotonics has been the subject of numerous recent theoretical and experimental 
works [4||6p31J . For example, two identical (symmetric) waveguides or membranes can attract 




Fig. 1. Schematic of single-membrane (asymmetric) structure: a photonic-crystal (holey) 
membrane (thickness hi = 0.2a) consisting of a square-lattice of air holes (radius R = 0.2a) 
on silicon is suspended (separation s) on top of an unpatterned (homogeneous) silicon slab 
(thickness h2 = 0.2a) sitting on top of a semi-infinite silica substrate. Light is incident on 
the membrane from the normal direction (top). 



or repel depending on the phase and frequency with which they are excited |[5][T4]|32j|35|, and 



this prediction led to several experimental realizations l|4]|I3[?9ll?3- Here, we extend this con- 
cept to greater generality by demonstrating that tunable attractive and repulsive forces can also 
arise in highly asymmetric structures coupled to external radiation, combining a suspended ho- 
ley membrane with an unpatterned layered substrate, forming two coupled planar waveguides as 
shown in Fig.[T] Such asymmetric geometries have a clear advantage in ease of fabrication and 
experimental characterization compared to symmetric double-membrane structures (in which 
two membranes must be both patterned and suspended). Although our system is asymmetri- 
cal, repulsive/attractive forces can still be thought of as arising from the bonding/anti-bonding 
level repulsion of guided-wave resonances similar to the phenomenon studied in symmetric 
systems 0^3|24j. In contrast to the symmetric case, however, it is more difficult to identify 
precisely which resonant modes of the large-separation (isolated) waveguides participate to pro- 
duce a resonant force, especially at higher frequencies, since the absence of exact degeneracies 
means that many more modes can potentially interact. Unlike symmetric systems, however, we 
can obtain bonding/antibonding forces that either increase or decrease with separation and more 
generally are nonmonotonic. We analyze the resonant forces from the perspective of perturba- 
tion theory. In unpatterned multilayer systems, we prove that only repulsive (not attractive) res- 
onances can occur. Although experimental characterization of forces in any membrane system 
requires accurate measurement of the membrane displacement, we point out that the transmis- 
sion/reflection resonances of the membrane itself can be used as a sensitive position sensor. 
Finally, in contrast to previous theoretical works that relied primarily on frequency-domain cal- 
culations, we propose a time-domain stress-tensor computational method that obtains the entire 
force spectrum from a single calculation via Fourier transforms of a short pulse. 

Optical forces can arise due to radiation pressure or gradient forces |3, 36|. Radiation pres- 
sure can be thought of as exchange of momentum between a photon (momentum hco/c) and 
matter, and as a consequence it is easily seen that light with incident power P exerts a force 
F = P/c on a. planar surface if the light is 100% absorbed or a force F = 2P/c if itis 100% 
reflected. Therefore, the ratio Fc/P is a. useful dimensionless measure of the strength of an 
optical force. Gradient forces, as shown recently |[5], can arise from the evanescent interaction 
between localized optical modes, and the resonant increase in the field intensity greatly en- 
hances the force for a given input power [31j, so that \F\c/P 1. Such large forces enable 
strong, tunable optomechanical interactions |13[ | T5}[22l [28| [36l [38l[39l , which have applica- 



tions such as optical cooling [ l|8[ |12[|4Q|^2 l, optical tweezers and traps |[3l [43|j47| , and optical 
switches |[9p7|. Successful demonstrations of other prominent resonant optomechanical effects 



include coherent mechanical oscillation (amplification) of mesoscopic objects with long vibra- 
tional lifetimes, and demonstrations of the photon-phonon strong-coupling regime via dressed 
optical states |[28l[g|48j|49l. 

We showed in Ref. ||5| that if two identical waveguides or resonant cavities are bought to- 
gether, the mutual interaction of these (degenerate) resonances or guided modes can induce 
a splitting of the modes into pairs characterized by attractive and repulsive mechanical forces, 
analogous to the well-known bonding/anti-bonding states formed by the level splitting (avoided 
crossings) of interacting degenerate states in quantum systems pQ| . In submicron-scale pho- 
tonic devices, these forces are strong enough to yield displacements and other mechanical ef- 
fects that have been observed experimentally ||4|p^[20)[26|. Furthermore, the frequency and/or 
phase of the optical excitation can be controlled to yield tunable optomechanical effects, even 
switching the sign of the force from attractive to repulsive. There are two ways to induce repul- 
sive and attractive interactions, depending on whether the incident power comes in the form of 
a guided mode or external radiation (the focus of this paper). First, one can inject light parallel 
to the membranes/waveguides, exciting guided modes that propagate along the waveguides and 
interact evanescently. In this case, the sign of the interaction is controlled by the relative phase 
of the modes in the two waveguides |[5]|T4||34[|35|. A similar effect occurs by controlling the 
relative phase of two coupled cavities (e.g. microspheres or microdisks), but in this case the 
bonding/anti-bonding resonances also have different frequencies that can be used to control the 
sign of the force ||20 37j 51 1. It is also possible to design asymmetric waveguide/cavity struc- 
tures (e.g. a dielectric waveguide and a microdisk resonator) with repulsive and attractive inter- 
actions as long as both structures support propagating modes, again with light incident along 
the waveguide direction |j4]|30|. (Coupled propagating modes in asymmetric geometries were 
also recently shown to lead to non-monotonic forces |52|.) Second, one can shine light per- 
pendicular to the membranes; if the membranes are perforated by periodic holes (or any other 
periodic modulation), normally incident radiation can couple via diffraction to guided-mode 
resonances within the membranes, which again couple evanescently. In this case, because the 
bonding/anti -bonding resonances have different frequencies, the sign of the force can be con- 
trolled by the frequency of the incident light. (By considering lateral shifts, one can also obtain 
lateral forces and other effects |24|.) As the periodic modulation (e.g. the hole radius) is made 
smaller, the lifetime (or quality factor Q) of the guided-wave resonances increases |24|, the 
resonant fields become stronger (intensity ~ Q), and thus the resonant forces become stronger 
~ Q (albeit narrower in bandwidth ~ 1/2). This force enhancement is ultimately limited only 
by losses (absorption or scattering from finite size or disorder). Another limitation is that the 
narrow bandwidth translates into a sensitive dependence of the force on the separation of the 
two membranes (since the resonant frequency shifts with separation). 

As indicated in Fig.[T] this paper considers the case of normal-incidence light on a suspended 
membrane, but in contrast to previous works, we only have a single suspended membrane (e.g. 
silicon) over a solid unpatterned substrate. The substrate must still support a guided mode of 
its own in order to obtain forces by evanescent coupling, and in our case this is achieved by a 
layered structure of a higher index layer on a lower-index substrate (e.g. silicon on silica). As 
in previous works, the membrane is periodically patterned to couple resonantly with normal- 
incidence light, and thanks to evanescent coupling this periodicity is also sufficient to couple 
light to guided resonances in the unpatterned substrate. Because the two waveguides in our 
case are so different, however, the interaction is more complicated than the degenerate level- 
splitting that occurs in symmetric systems, leading to nonmonotonic force dependences as well 
as transitions in the sign of the force. The separation and frequency dependence of the force in 
membrane structures has been previously exploited for achieving a variety of optomechanical 
effects. For example, one can obtain mechanical oscillators with dynamically tunable "spring 



constants," even flipping the sign to yield an unstable equilibrium and mechanical bistability, 
with a tunable relative strength of the linear and nonlinear terms (with arbitrarily strong relative 
nonlinearity possible if the linear term is canceled) 1 13 p5p2p8p6l[38p9| . Here, we point out 
that by operating with light consisting of two frequencies, rather than a single frequency, one 
can tailor the spring constant of the membrane without changing its mechanical equilibrium 
separation. Furthermore, we argue that, with appropriate design, it should be possible for repul- 
sion to dynamically activate at small separations, creating a feedback mechanism for combating 
stiction arising from other forces (e.g. electrostatic ||53j[54| or quantum ||55j|56| interactions). 

2. Computational Method 

Previous numerical methods for computing optical forces have mainly focused on frequency- 
domain approaches, with some exceptions |57}j59| . Given two objects separated by a distance 
J, the force on one of the objects can be computed in one of at least two ways: One approach 
involves computing the derivative of the eigenfrequencies CO of two membranes separated by 
distance ^ as a function of s, which can be related to the force via the relation F = l/co dco /ds^ 
[37j, where F >0 corresponds to repulsion. Ref. |60| generalized this formula to handle the case 
of resonances with finite lifetimes 2n(0, rigorously justifying our earlier suggestion [ 37 1 that 
the change in lifetime with separation has a negligible (higher-order) contribution to the force 
on resonance. This approach is problematic, however, in general circumstances where there 
may not be well-defined resonant modes with negligible loss, nor does it include cases where 
there is a superposition of the resonant mode with other waves (e.g. light from an external 
source). Another approach involves computing the force via an integral §^(T) • dS of the time- 
average Maxwell stress tensor 



(1) 



around some bounding surface S lying in vacuum |2|. For resonant modes with negligible loss, 
{Tij) can be computed directly from an eigenmode calculation [23 1. More generally, including 
the case where the fields are excited from an external source whose effect must be included 
from some time-harmonic current source i{x)e~^^\ one can solve a set of linear equations for 
the resulting time-harmonic fields E and H by a variety of methods (e.g. finite elements or 
differences in the frequency domain, or transfer-matrix methods) |61, appendix D], and use 
these fields to compute Tij | [T0|[26[|62 -64|. This approach, however, has the drawback that if the 



force at many frequencies (0 is desired, one must perform many separate calculations (one for 
each (O). 

If a broad-band force spectrum is desired, an attractive alternative is to compute the stress 
tensor via the Fourier transform of a short pulse in the time domain (e.g. finite-difference time- 
domain, FDTD 1 65 1), yielding the entire frequency spectrum at once. Here, one simply evolves 
Maxwell's equations in response to a pulse source [e.g. ~ J(x)5(^)] in time, accumulating the 
discrete-time Fourier-transform [f{(o) ~ £„/(^A^)^^^^A^] of both the electric E and magnetic 
H fields over the stress-tensor surface S, and at all desired frequencies (0 | 66 67 1. These Fourier- 
transformed fields then yield the stress tensor and hence the force. Of course, the force must be 
normalized in some way, and here the dimensionless Fc/P normalization is very convenient. 
One simply does a separate calculation, with no structure (vacuum), to compute the Fourier- 
transformed incident fields and hence the incident power P{co). Dividing F{co)c/P{(o) yields 
the dimensionless force spectrum, where all arbitrary normalization factors (e.g. the incident 
pulse spectrum or the normalization of the Fourier transform) have canceled. (Matters are more 
complicated in a nonlinear system, of course.) 




Fig. 2. {Left:) Resonance frequency ft) (units of Inc/a) a function of separation s (units of 
membrane period a), for both the single-membrane (asymmetric) structure of Fig.[T](soUd 
Unes) as well as the double-membrane (symmetric) structure of Ref. ||24| (dashed lines). 
The insets show the electric field component Ex in the x-z plane {y = 0) near ft) at a particu- 
lar s = 03a. In the symmetric case, the attractive and repulsive modes are in-phase and out 
of phase, respectively, as expected. (Right:) Resonant (peak) force Fc/P (units of incident 
power P/c), at the resonant frequencies ft) plotted on the left figure, as a function of s. The 
bottom inset shows the broad-bandwidth force spectrum of the asymmetric structure at a 
particular s = 0.2a, showing both the bonding (F > 0) and antibonding (F < 0) resonances. 
The inset also denotes what is meant by resonance frequency ft) and peak force F. 



In what follows, we exploit our FDTD approach to compute forces on the geometry of Fig. [T] 
All of the subsequent calculations were performed using MEEP, a free FDTD simulation soft- 
ware package develped at MIT |67|. We find that discretization errors coming from our finite 
resolution of 40 pixels/a affect the computed force spectra by no more than a few percent. 

3. Membrane Forces 

In this section, we explore attractive and repulsive resonances in the asymmetric membrane 
structure of Fig. [T] along with possible applications and the underlying theory. 

3.1. Symmetric and asymmetric systems 

We begin by considering two related membrane structures: the asymmetric structure of Fig.[T] 
with a perforated silicon membrane over a silica substrate; and also a symmetrized version 
consisting of two identical perforated silicon membranes, examined previously by Ref. |24|. In 
both cases, the membranes are illuminated from above by normal-incident {z direction) light 
polarized in the x direction, and we consider the resonant frequencies and the resulting forces 
as a function of frequency and separation. 

In the symmetric case, it is well known that each resonance of the individual membranes 
splits into two resonances of the coupled two-membrane system: "bonding" and "anti-bonding" 
modes, in which the individual resonances are excited in phase and out of phase, respec- 
tively |5|. The frequencies of these two resonances as a function of membrane separation 
are shown as dashed lines in the left plot of Fig. [2j and as expected the frequency splitting 
vanishes as the separation increases (and hence the membrane coupling decreases) |5|. The 
corresponding field patterns are shown as the right two insets of the left plot, and display 
the expected phases. Each resonant mode corresponds to a resonant peak in the optical force, 
and this peak force (at the resonant frequency) is plotted as a function of separation as dashed 



lines in the right part of Fig. |2] As expected, the bonding and anti-bonding modes correspond 
to opposite-sign attractive and repulsive forces between the membranes, respectively, and the 
force becomes stronger as the separation decreases (increasing the membrane interactions). The 
attractive force in the bonding case has slightly larger magnitude than the anti-bonding repul- 
sion, which can be explained by the larger field overlap in the former case due to the lack of a 
node in between the membranes | [33l[37| . 

In the asymmetric case, the mode of the isolated membrane is not the same frequency as the 
corresponding guided mode of the isolated layered- substrate structure (although the parameters 
can be adjusted to force a degeneracy if desired). The mode of the silicon on silica (oxide) sys- 
tem is actually a lossless waveguide mode (lifetime ^ Q = o^), not a resonance; it is only when 
the membrane is brought into proximity with the oxide that the membrane's periodicity a allows 
guided modes at wavevector 2n/a (and multiples thereof) to couple to normal-incident radia- 
tion | [6T| . In this case, the layered- substrate (waveguide) mode that is nearest in frequency to 
the isolated-membrane resonance frequency is the lowest-order waveguide mode of wavevector 
In /a. Because the two mode frequencies are no longer degenerate, when the resonant frequen- 
cies of the asymmetric case are plotted versus separation as solid lines in the left part of Fig.[2j 
the frequency splitting no longer vanishes as the separation increases. Nevertheless, there is a 
frequency splitting or "level repulsion," explained below in terms of second-order perturbation 
theory, which becomes significant for small separations, and the corresponding field patterns 
display the qualitative phase characteristics of bonding/anti-bonding modes (insets). As a con- 
sequence, as considered theoretically below, the force spectrum in the asymmetric case (right, 
inset) indeed displays the characteristic attractive and repulsive resonant peaks of bonding/anti- 
bonding modes. The peak force (at resonance) versus separation is plotted as solid lines in the 
right part of Fig. [2j and has similar sign as in the symmetric case. Of course, the system is 
now more complicated than the symmetric case in a variety of ways (e.g. the field patterns are 
no longer symmetrical/anti- symmetrical and the lifetimes as well as the frequencies depend 
strongly on separation), so the peak force versus separation dependence is significantly differ- 
ent: First, the peak bonding (attractive) force decreases as s decreases and reaches a constant 
value as ^ ^ oo. Second, the ratio of the antibonding (repulsive) to bonding force becomes in- 
creasingly larger at smaller separations (e.g. it is more than a factor of 2 larger at ^ = 0.1a), in 
contrast to what is normally observed | 33 , 37 1. 

To understand these features of the force, we use the fact (reviewed in Sec. |3.3| ) that the force 
is proportional to both the separation {s) dependence d(o/ds of the frequency and also the life- 
time Q. As discussed below, perturbation theory indicates that d(o/ds of a nondegenerate mode 
in one object decreases proportional to the square of its field overlap with the other object | 68 1. 
For a mode which as ^ ^ oo approaches a leaky mode of the perforated membrane — in our 
case, the antibonding mode — the lifetime Q asymptotes to a nonzero constant and hence the 
product Qd(o/ds 0; correspondingly, the force tends exponentially to zero with s. Similarly, 
the force must tend to zero for both the bonding and antibonding modes of a symmetric mem- 
brane (where all modes are leaky as ^ ^ oo). On the other hand, for a mode that asymptotes 
as 5- ^ oo to a lossless guided mode of the unpatterned substrate — in this case, the bonding 
mode — the lifetime Q diverges sls s ^ oo. In fact, perturbative scattering theory |68 | indicates 
that the scattered power, and hence 1/Q |[6T|, is proportional to the square of the field overlap 
with the periodic membrane (the source of the scattering loss). Hence Q diverges at the same 
rate at which d(o/ds vanishes, and thus the force should asymptote to a nonzero constant as 
^ ^ oo. These behaviors are precisely what is observed in Fig. |2] (right): the peak forces of the 
symmetric system and the asymmetric antibonding mode decrease monotonically to zero with 
increasing s, while the peak force of the asymmetric bonding mode increases monotonically to 
a constant with increasing s. The corresponding variation in Q is shown in Fig. [3] (left). Note 
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Fig. 3. {Left:) Optical force Fc/P on the single-membrane structure of Fig. [T] as a function 
of the frequency CO of incident light of power P, for various separations s. The insets show 
typical Ex field patterns (in the x-z plane, at j = 0) for both the attractive (left) and repulsive 
(right) resonances. (Right:) Optical force Fc/P as a function of separation s for incident 
light input at various frequencies CO G [0.48,0.5] (Inc/a). The bottom inset shows Fc/P 
for light input over a lower frequency range CO G [0.41,0.43] Inc/a. The force versus s 
plot was obtained by fitting the force spectrum obtained via FDTD at a few 5 to a sum of 
Lorentzian resonances, and then interpolating the resulting Lorentzian parameters over a 
denser s range. 



that for s < 0.2a, the Q of the antibonding mode increases rapidly with decreasing s, leading 
to an increasing antibonding force that is many times larger than the corresponding bonding 
force. In a practical system, the behavior of the bonding mode will be further modified by the 
presence of loss in the isolated- substrate guided mode (from finite-size effects, roughness, ab- 
sorption, etcetera) — this will cause its Q to saturate to a finite value. In this case, the force will 
behave nonmonotonically. it will initially increase, but will then decrease to zero as s goes be- 
yond the saturation point of Q (while dco/ds continues to decrease). Thus, the lifetimes of both 
the membrane and substrate could be exploited to tailor the s dependence of the force in this 
and other similar systems. 

3.2. Tunable mechanical properties 

Because the frequency and magnitude of the resonant forces change as s is varied, it is interest- 
ing to study also the separation dependence of the force for light incident at a single frequency, 
which alters the mechanical dynamics of the membrane. Here, we consider the effect of light 
incident at a single frequency, and then extend our analysis to the case of two frequencies. 
(More generally, it may prove interesting to study the dynamics of the membrane for modu- 
lated pulses.) Figure [3] plots the optical force Fc/P as a function of separation s for incident 
light at various frequencies (0. As expected, the lifetime ~ g of the attractive peak, with 
concentrated in the layered substrate (shown on the inset), becomes infinite (2 ^ 00) as ^ ^ 00 
due to the reduced coupling between the infinite- 2 substrate mode and the finite- 2 PhC reso- 
nance. At a fixed frequency, changing s can move the system into or out of resonance, leading 
to a dramatic ^-dependence of the force. (Indeed, the s dependence can be much sharper than 
shown here, e.g. if the hole diameter is shrunk to increase the Q of the resonances.) One can 
obtain transitions in the sign of the force, from attractive to repulsive and vice versa, as s is var- 
ied, leading to stable and unstable equilibria, not yet including the mechanical restoring force 
from the membrane supports, and even multiple equilibria. 



When mechanical forces are included, two things can happen. If the optical force is nonzero, 
the mechanical equilibrium point of the membrane will shift and the slope of the force curve (the 
spring constant Ko = dF/ds) will be altered. If one operates at a point where the optical force 
is zero, then the equilibrium position is unaltered but the spring constant is changed. The total 
spring constant, including linear mechanical restoring forces on the membrane, will be given by 
K= Ko-\-Kjn, where denotes the mechanical spring constant. Whereas is frequency- and 
power-independent, Ko exhibits a very sensitive dependence on both, and therefore by choosing 
CO and the incident power it is possible to tune the total spring constant of the system | [69}j72| . 
On the one hand, if one chooses CO so that Ko/Km> 0, then optical forces act to increase K and 
therefore stiffen the stable equilibrium. In systems driven by undesirable thermal fluctuations, 
this effect has been exploited for "cooling" the resulting vibrations |^0|^0j|^. On the other 
hand, if one chooses CO so that Ko/Km < 0, then K can be decreased and even flip sign as the 
optical power increases, leading to an unstable equilibrium and bistable behavior |69j[7T|. Near 
an exact cancellation /c^ ^ — /c^, the linear term in the ^--dependence of the force is decreased 
relative to the higher-order nonlinear terms (which include both optical and mechanical terms), 
allowing arbitrarily strong nonlinear mechanical effects, and even a strictly nonlinear regime 
of operation (Ko = —k^) where effects like bistability, hysteresis, and frequency conversion 
should be readily observable |36, 74|. 

For light incident at a single frequency CO, sign transitions in the force occur as the structure 
moves toward or past a force resonance, as illustrated in Fig. |3] Away from these resonances, 
the "background" dimensionless force Fc/P is attractive and bounded from above by 2 (see 
Sec.|4]), leading to negligible optical spring constants (small Ko) at the corresponding equilibria 
separations [see Fig. [3] (inset)]. This generally does not preclude a strong modification of the 
mechanical properties of the membrane (achieving large Ko ^ K^) since it is also possible to 
operate at separations where Fc/P is large and has linear slope (dF/ds ~ s), although this in- 
evitably causes a change in the initial mechanical equilibrium separation of the membrane |39|. 
For applications in which achieving a large Ko without modifying the initial mechanical sepa- 
ration (i.e. achieving a large dF/ds at a position where Fc/P = 0) is important, then a different 
scheme is required. For example, rather than operating with incident light at a single frequency 
CO, one can instead consider the combined effect of light incident at two different frequencies 
(0+ and CO-, near the attractive (bonding) and repulsive (antibonding) resonances, respectively. 
This idea is illustrated in Fig. |4] (left), which shows the optical force Fc/P as a function of sep- 
aration s for incident light of power P = P^-\-P- consisting of two frequencies, [chosen in 
the region (0+ G [0.41,0.424] {27rc/a)] and CO- = 0.495 i/lnc/a), of corresponding power P+ 
and P-, respectively. From Fig. [3] it is clear that incident light at CO- leads to a repulsive peak 
dii s ^ 0.2a whereas incident light in the (0+ range leads to attractive resonances in the range 
s e [0.1, 0.3] a. As the attractive and repulsive peaks excited by CO- and co^ come close to one 
another, the transitions in the sign of the force become more pronounced, leading to larger Ko. 
To quantify the enhancement in the spring constant. Fig. |4] (right) plots the absolute value of 
the optical spring constant \ Ko\ (units of P/ca) as a function of 0)+, for different values of the 
ratio 7] = P-/P^ of power in CO- versus (0+, where dashed/solid lines correspond to negative 
(unstable) and positive (stable) Ko, demonstrating orders -of-magnitude enhancement in Ko. For 
example, the peak spring constant \Ko\ in the case where r] = 1 is ^ 10"^ {P/ca), whereas it 
is smaller than 10 (P/ca) in the case of incident light only at 00- , corresponding to r] ^ oo. 
An alternative scheme that allows tailoring of the optical spring constant near equilibrium was 
explored in Ref. | [5T| , although in that case the effect is achieved by the presence of multi- 
ple bonding/antibonding pairs in which opposite-sign force resonances were designed to occur 
at closely spaced frequencies, whereas here there is no need for the resonances to be closely 
spaced. 




Fig. 4. {Left:) Optical force Fc/P as a function of separation 5, for light incident at two 
frequencies ft)+ (varied) and ft)_ = 0.495 {Inc/a), with corresponding power P+ and P_, 
respectively. Dashed lines show the force for P_ = 0. (Right:) Absolute value of optical 
spring constant \ Ko\ (units of P/ca) as a function of frequency ft)+. Dashed and solid lines 
correspond to negative (unstable) and positive (stable) values of /Q,, plotted for different 
values of rf = P-/P^. Both Fc/P and Kq are normalized against the total input power 



Additional forces on the membrane arise at small separations due to residual static charges 
and also due to quantum/thermal fluctuations (Casimir and van der Waals forces), which are typ- 
ically attractive |75 -78 1 and may lead to "stiction" problems in micromechanical (MEMS) sys- 
tems where moving parts are forced into contact | [55l56|79j | . Here, the separation dependence of 
the classical optical force can potentially be used to combat such stiction effects |80|8T1. Not 



only can one exploit a repulsive resonance to oppose stiction, but the separation dependence 
means that such a repulsion can be designed to take effect only if s inadvertently falls below 
some threshold. That is, a repulsive resonance for small s can be used as a feedback effect to 
reduce the chance of stiction without significantly altering the mechanical dynamics at larger s 
where the incident light is out of resonance. In an upcoming manuscript, we will demonstrate 
how these effects can also be exploited to design integrated, all-optical and accurate tech- 
niques for measuring the Casimir effect that rely on measuring static displacements rather than 
forces or force gradients. 

Figures|3]and|4]only show a small sample of the kinds of optical effects that can be observed 
in evanescently-coupled systems. In particular, there are many degrees of freedom and possi- 
bilities to explore, especially if one is not restricted to symmetric structures or operating at a 
single frequency. For example, the magnitude of the optical forces shown in Fig. [3] are by no 
means the largest possible, since larger forces can be obtained merely by increasing Q at the 
expense of bandwidth (and s insensitivity). Even more complicated behaviors can be obtained 
by increasing the number of resonances, and the choice of resonance offers a corresponding 
choice of lengthscales or operating frequencies. 

3.3. Level repulsion in asymmetric membranes 

For well-defined (long lifetime) resonant modes, perturbation theory has been used to analyze 
the relationship between the force and the resonant frequency/lifetime |[34 37j72|, and this 



relationship can also be used to illuminate the relationship between the sign of the force and 



^The interaction of normal-incident light with the membrane in this system can be exploited to simultaneously 
control and measure the membrane's equilibrium separation. 



the field distribution, as well as the physical origin of resonant repulsion. 

As long as the interactions in the system are dominated by a discrete set of well defined leaky 
resonances (lifetime 1 /period), then one can apply standard methods of discrete- spectrum 
time-independent perturbation theory |68|. (For non-resonant systems with a continuum of 
non-localized modes participating at every frequency, perturbative methods are much more 
complicated |82|.) In this case, we proceed in three steps: first, we relate the force F{s) to 
the ^-dependence dco/ds of the resonant frequency CO (generalizing a previous result |37|); 
second, we connect dco/ds to the electric-field distribution of the resonance; third, we incor- 
porate the level-repulsion effects of nearby resonances via second-order perturbation theory. 
There is one important complication that arises in resonant systems which does not arise in the 
lossless guided modes where dco/ds effects were previously derived: the quality factor Q (life- 
time Q/ CO) of the resonances may in general also vary with separation. Ref. | 60 | considered 
the effect on the force due to changes in both CO and Q and found generally that at resonance 
(including the structure considered here), dQ/ds effects have little if any effect on the force (be- 
ing higher order in 1/g, as we previously suggested without proof |^37J). We therefore neglect 
dQ/ds in the following, in which case the expression for the force is: 

Fc ^ _ c ^ f Q\ ^ Q^dco 

P Pds \ co) co^ ds' ^ ^ 

where F <0 denotes a repulsive force. Although Eq. ^ is derived rigorously from coupled- 
mode theory in Ref. | [60| , a simple justification for the same result can be obtained if one 
neglects radiative or absorptive losses, to treat the resonator as a closed system (i.e. equivalent to 
changing the separation slowly compared to 1 /(O but quickly compared to the lifetime). Given 
an incident power P, the energy U stored in a resonance of real frequency CO and quality factor 
Q is (by definition of Q | 61 1) given hyU = PQ/ CO. Neglecting radiation loss, any change dU /ds 
in the energy must be due to a mechanical force F = —dU /ds, resulting in Eq. Equivalently, 
F /U = 1/(0 {dco/ds), a result we previously derived for guided (Q = oo) modes |5 1. 

The dependence of CO on s can be predicted by perturbation theory. In particular, the first- 
order change 5o)^^^ to the frequency CO coming from a small change Ae in the permittivity of a 
system with original permittivity e is readily expressed as: 

5coW _ l(Ea,|Ae|Ea,) 



CO 2 {Eco\£\E 



CO/ 



In this case, however, Ae is not small: at a given point near the interface, £ is changing dis- 
continuously as that interface moves past the point. In this case, perturbation theory must be 
derived more carefully | [68l . For an interface from £i to £2 that is moving by A.s' (towards £2), 
assuming isotropic materials, the numerator of Eq. ^ changes to: 

(E|Ae|E)^(E|||A.(ei-e2)|E||) 

-{D^\As{£^'-£^')\D^) (4) 

Without loss of generality, we can hold the upper membrane fixed and move the substrate (or 
lower membrane) away by A^. From Eq. (|4]), the way to obtain attractive (dco/ds > 0) and 
repulsive (dco/ds < 0) resonant effects is clear. If we hold one membrane fixed and move the 
substrate (or the other membrane), dco/ds will be positive (attractive) if |Ep is peaked where 
£1 < £2, i-e. on the air/silicon interface (adjacent to the upper membrane). Conversely, dco/ds 
will be negative (repulsive) if |Ep is peaked where £1 > £2, i.e. on the silicon/oxide interface 
(away from the upper membrane). Precisely such field patterns can be observed in the insets of 
Fig. |2] and Fig.|3j the repulsive and attractive modes have Ej peaked at the expected interfaces. 



Note that a homogeneous substrate, e.g. semi-infinite silicon or oxide, has no interface except 
for the air interface adjacent to the upper membrane, so in this case a repulsive force cannot 
arise by this mechanism, as demonstrated for the h2 = case on the inset of Fig. [5] below. An 
exception to this rule is discussed in Sec. [4j in which repulsive forces arising from radiative 
modes are analyzed, a situation where a lack of normalizability causes the perturbation theory 
to break down. 

The above discussion indicates which field patterns might be expected to lead to repulsion 
and attraction, but does not explain how such field patterns can arise. For the case of a sym- 
metric membrane, symmetry considerations predict that the degenerate modes of two isolated 
membranes will split into even/odd bonding/anti-bonding pairs by degenerate first-order pertur- 
bation theory [ [83| . The presence of a nodal plane bisecting the anti-bonding mode (assuming 
its field is dominated by Ejy and not E^) means that the field pattern will be stronger on the 
far sides of the membranes, leading to a repulsive interaction as observed in Fig. [2] and Fig. [3] 
and as predicted above. For asymmetric membranes, however, there is typically no degener- 
acy, and the corresponding two-mode interaction must instead be analyzed by second-order 
perturbation theory, which plays a role for sufficiently small separations (large interactions). 
When two isolated waveguides each have a mode with nearby frequencies, and one brings the 
waveguides together so that the mode fields overlap, second-order perturbation theory predicts 
a contribution to Aco that tends to split the two frequencies: 

gcQp) ^ 1 |(Ea,|A£|Ea,)p 1 ^ / 0)^ \ | (E^. | [Eg,) ^ 
CO 4 \{E^\e\E^)\^ 2j^^\(o'^-co^) (E^,|^|E^,) (E«|£|Ea,) 

(As above, the overlap integrals are modified into Ey and components for motion of dis- 
continuous interfaces |[68|.) Note that the (o'^ — co^ term pushes CO away from co\ and be- 
comes stronger as the frequencies become closer (with the most dramatic case being degener- 
ate modes, where the derivation is modified). Because of the competition between the first- and 
second-order terms, which may be comparable in magnitude for small separations (large over- 
laps), it is possible for the force near a resonance to switch signs with separation, a possibility 
that is demonstrated in Sec. l3.4l below. 

3.4. Multi-modal interactions 

Previously, we considered resonances at a relatively low frequency (compared to 27lc/a), where 
the only relevant interactions were between two modes (one for each isolated membrane or 
substrate). At higher frequencies, the density of states generally increases, and thus many more 
resonant modes are typically present. Correspondingly, the inter-modal interactions become 
more complicated, and it is not always possible to identify individual pairs of bonding/anti- 
bonding modes. However, the qualitative features of repulsive and attractive resonances are 
still present, although the additional complexity provides more degrees of freedom leading to 
more complicated force phenomena. 

For example, the force and reflection spectra for the asymmetric membrane system of Fig. [T] 
are shown in Fig. [5] in a higher frequency window (about double the frequencies in Fig. [3]). As 
before, the force spectrum (left) shows both repulsive and attractive resonances. In this case, 
however, we actually observe a force resonance changing sign as a function of separation, which 
physically can be interpreted as different terms dominating in the perturbation theory [Eq. ([5])] 
at different separations. In the reflection spectrum (right), these resonances correspond to Fano 
shapes (adjacent peaks and dips), a well-known consequence of the coherent combination of a 
resonant process with direct transmission through the slabs ||9[|24[|84| . Because the reflection 
spectrum depends sensitively on the separation, the peak locations from a broad-bandwidth low- 
intensity source could be used to accurately determine the separation in an experiment. Note 
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Fig. 5. {Left:) Optical force Fc/P as a function of frequency (O for light of power P incident 
on the single-membrane structure of Fig. [T] for various separations s. The bottom inset 
shows the force (solid lines) and reflection (dashed line) of the same geometry but for 
h2 = and s = 03a. (Right:) Corresponding reflection spectrum P as a function of ft). The 
open circles indicate frequencies for which there exist force minima or maxima. The insets 
show the electric field component Ex in the x-z plane (y = 0) at a particular s = 0.2a, and 
at the indicated frequency points CO = 0.64{2nc/a) (left) and co = 0.681 (27rc/a) (right). 

also that the resonant peaks in this frequency window approach one another as the separation 
decreases, which means that the largest contribution to level repulsion in this case is coming 
from interactions with other modes (outside this frequency window, not shown); this is verified 
by examining the field patterns (upper insets), which clearly correspond to completely different 
modes in the membrane and not just a relative-phase change. As in Sec. |3.2[ fixing the frequency 
and plotting the force versus separation reveals a force that changes both magnitude and sign 
as a function of separation. 

In the absence of a guided mode in the substrate, the physics of this situation is completely 
changed — one no longer has level repulsion effects, since there are only resonances in the 
membrane. A general argument for the attractive nature of the resonant force is given in the 
previous section. This occurs, for example, if the oxide substrate is replaced simply by a low- 
index (oxide, n = 1.5) substrate, which supports no guided modes of any kind on its own. 
However, the resonant frequencies (and lifetimes) of the membrane are still modified by the 
proximity of the substrate, so there is still a resonant force, and in this particular geometry we 
find that the resonant force is always attractive. Such attractive forces can however be exploited 
in integrated photonic devices as a way of tuning the mechanical response of the devices | 52j 
[85| . The force and reflection spectra for this case are plotted in the middle inset of Fig. [5] (left), 
for a particular separation s = 03a. 

4. Fabry-Perot Forces 

Our focus thus far has been on periodic structures supporting resonances whose coupling leads 
to both enhanced attractive and repulsive forces. However, at large separations compared to 
the evanescent tail of the participating guided modes, evanescent bonding/anti-bonding effects 
lead to negligible force enhancement. On the other hand, there exists an alternative and comple- 
mentary force-enhancement mechanism that can play a role at large separations: light normally 
incident on two separated planar objects can be resonantly enhanced due to the Fabry-Perot 
cavity formed between the objects by reflections from adjacent surfaces of the objects, and the 
extent of this enhancement will increase with the reflectivity of the objects (i.e. with the Q of the 



Fig. 6. Schematic of system consisting of two multilayer objects [labeled as (1) and (2)] 
separated by a distance s. A two-dimensional cross-section for the particular case of two 
quarter- wave stack mirrors with a defect (yellow) is shown on the right. 



cavity). Although such Fabry-Perot force-enhancement mechanisms have been considered in 
previous work [20^74^86 - 89] , in applications ranging from gravitational-wave detection | 90 | to 
optical cooling | |42||91f|94| , it is interesting to explicitly compare the two resonant mechanisms. 
In this section, we briefly consider the kinds of repulsive and attractive forces that can arise in 
systems consisting of unpatterned multilayer objects, emphasizing some of their similarities 
and differences compared to forces arising from evanescently coupled resonances. 

Resonant radiation pressure within Fabry-Perot cavities has a long history, dating back to 
work in the 1960s on interferometer sensitivity [95J , and has since been considered both the- 
oretically and experimentally for many applications |[36j[88|, such as nonreciprocal phenom- 
ena 1 96 1, optical cooling | [92}j94| , and tunable optical springs f47|. If the space between two 



partially reflecting mirrors (e.g. Bragg mirrors) is viewed as a waveguide, then the resonance 
frequency for normal-incident light corresponds to a slow-light (zero group- velocity) band edge 
where radiation pressure is enhanced fTO",^!. In all of these cases, the pressure is repulsive, 
as one might expect for light bouncing between the two objects (and was argued in general for 
two semi-infinite objects 1 64 1). We show this in general below, for radiating modes (not guided 
modes) and any unpatterned multi-layered structure (in the absence of gain). 

Here, we consider a general class of geometries, depicted in Fig. [6j consisting of two 
unpatterned (translationally invariant in two directions) planar multilayer objects separated 
by distance s in vacuum, denoted as objects (1) and (2), characterized by complex reflec- 
tion/transmission amplitudes r\/t\ and r2/^2 at a given frequency, respectively (satisfying 
ky^P + ky^P = 1 in the absence of absorption or gain). The dimensionless force Fc/P on ob- 
ject (1) due to light normally incident from above at frequency CO can be readily computed 
(using a simple transfer-matrix analysis to obtain the stress tensor [97 J ) to be: 



^ = 1 

P 



ri^hr2e^'^F^ 



where =t\/[l — rir2exp(2/5)] is the induced field at the lower interface of object (1), and 
5 = Ittcos is the phase associated with the air gap. 

Elementary manipulation of Eq. ^ shows that the force is bounded above by Fc/P < 2, 
which is a bound on the attractive (positive) force. This result is a physical consequence of 
conservation of momentum: the light trapped between the objects can only act to repel them, 
whereas an attractive force can only arise from reflections from the uppermost surface. (The 
key difference compared to patterned membranes or guided modes is that the fields in between 
the objects are now purely propagating waves in which the field amplitude is proportional to 



the wave momentum |2|, with no evanescent component where these two quantities can be 
decoupled.) The maximum reflectivity is 100%, corresponding to Fc/P = 2. A related situation 
is one in which the lower object plays no role because its reflectivity approaches zero. In this 
case, for r2 = and 1^21 = 1, one recovers a standard result for the force on a single object, 
Fc/P = 1 + |ri p — 1^1 P ||95|, which is always positive (attractive) and is bounded above by 2. 
(Note, however, that all of these limits only apply in the absence of gain, which can alter the 
force by emitting additional photons from the objects |[98|.) On the other hand, the repulsive 
forces are unbounded, becoming arbitrarily large as | ri | and | r2 1 approach unity. This is simply a 
consequence of the repulsive pressure from the Fabry-Perot mode trapped between the objects, 
whose lifetime (and energy density) diverge in this limit. This familiar result has been exploited 
in many cavity-enhanced optomechanical systems as noted above. 

One can construct multilayer objects supporting exponentially localized resonances which 
couple to normally incident radiation. For example, this is the case if each object consists of a 
multilayer Bragg mirror with an embedded defect layer | |6T| . In this case, degenerate perturba- 
tion theory implies that two symmetric objects [such as those in Fig. |6] (right)], each with an 
identical embedded defect/resonance, should couple to form bonding/anti-bonding states where 
the resonances are in/out-of phase. Naively, one might suppose that this will lead to repulsive 
and attractive resonances, as in the coupled guided-mode case, but Eq. ^ indicates that this is 
impossible. The explanation is straightforward: although such defect modes are exponentially 
localized within the Bragg mirrors composing each object, they are propagating in the region 
between the objects where there are no mirrors (because the input beam is propagating in free 
space and no diffraction occurs). This means that the frequency splitting does not depend ex- 
ponentially on the separation between the two objects, and hence there is no resonant force 
enhancement via these modes, by the analysis of Sec. [3] 

5. Conclusion 

Optomechanical interactions are a rich subject of current research, and the use of evanescently 
coupled guided resonances enables an especially rich set of phenomena because of the pres- 
ence of both attractive and repulsive resonances. The ability to tailor and exploit guided reso- 
nances coupled via periodic modulations offers an exciting opportunity to procure complicated 
force effects at small separations. In this paper, we showed that functionality similar to that of 
previously studied symmetric-membrane systems can be obtained in asymmetrical membrane- 
substrate structures. From an experimental point of view, such asymmetrical structures are at- 
tractive in that only a single membrane need be suspended and patterned. From a theoretical 
viewpoint, because the resonant modes of asymmetrical structures need not come in degenerate 
pairs (unless degeneracies are forced), more than one pair of modes can have strong interac- 
tions, leading to the possibility of richer force phenomena |51 1. Correspondingly, the distance 
dependence of a force spectrum with multiple attractive and repulsive resonances can exhibit 
richer "optical spring" phenomena than are possible with repulsive resonances alone (as in 
Fabry-Perot resonances between mirrors); for example, one can operate at a zero of the optical 
force to tune the optical spring constant (in either direction) without altering the equilibrium 
position. 
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